Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Extended Star Formation in Dwarf Spheroidal Galaxies: the 
Cases of Draco, Sextans, and Ursa Minor 

■ C. Ikuta 1 '* and N. Arimoto 2 ' 3 
O ' 

o : 

£NJ ■ 1 School of Physics & Astronomy, University of Nottingham, Nottingham, NG7 2RD, U.K. 

' e-mail: chisato.ikuta@nottingham.ac.uk 

Q_l| 2 National Astronomical Observatory of Japan, 2-21-1, Osawa, Mitaka, Tokyo 181-8588, JAPAN 

■ e-mail: arimoto@optik.mtk.nao.ac.jp 

3 Institute of Astronomy, University of Tokyo, 2-21-1, Osawa, Mitaka, Tokyo 181-0015, JAPAN 

<n : 

| Received 18 October 2001; accepted 19 April 2002 

> ; 

Abstract. Star formation and chemical enrichment histories of the dwarf spheroidal galaxies (dSphs) Draco, 
Sextans, and Ursa Minor are investigated by means of chemical evolution models and a simulation code for 
colour- magnitude diagrams (CMDs) . The CMD simulation code is designed to fully consider effects of the chemical 
evolution on stellar evolution and photometric properties. For this aim, star formation and chemical enrichment 
histories are calculated consistently in the code. Comparisons between the chemical evolution models and the 
observed abundance patterns reveal that the star formation rates were very low (1 — 5% of that of the solar 
neighbourhood disc) and that the initial star formation continued for a long duration (> 3.9 — 6.5 Gyr) in these 
' dSphs. This star formation history can reproduce morphologies of the observed CMDs, such as narrow red giant 

Oh' branches and red horizontal branches and succeeds in solving the second parameter problem of the dSph Draco. 

Hence, both of the abundance patterns and the morphologies of the CMDs can be explained by the star formation 
5_l ■ histories characterised by the low star formation rate and the long duration of the star formation period. Because 

of the low star formation rates, plenty of gas remains at the final epoch of star formation. We suggest that gas 
C$ • stripping by the Galaxy results in termination of star formation in the dSphs. 
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1. Introduction various elements in the three dSphs Draco, Sextans, and 

Ursa Minor, which are all satellite galaxies of the Milky- 
Satellite dwarf spheroidal galaxies of the Milky Way have Way . They report that the abundance ratios of a elements 
been considered to be similar to globular clusters (e.g., (Q) gi) M& etc ) relative to iron ([ a / Fe ]) are nearly equal 
Hodge 1971), since they are composed of old metal-poor to golar at low metallicity ([Fe/H] < -1.4). This suggests 
stars and contain little gas. Recent analyses of the colour- that thc starg were formed in the gag en riched by Type 
magnitude diagrams (CMDs), however, reveal a surprising Ia supe rnovae (SNela) as well as by Type II supernovae 
fact that dwarf spheroidal galaxies (dSphs) are of compos- (SNe ^ that igj star formation m the dSphs continued 
ite stellar populations. Some of dSphs formed stars contin- Qver the lifetime (typically 1 - 2 Gyr) of SNe Ia progeni- 
uously and/or intermittently for many Gyr (e.g., the dSph torg Thig regult allows ug to direc tly study the chemical 
Carina: Smecker-Hane et al. 1996; the dSph Leo I: Gallart enrichment histories (CEHs) of the dSphs Draco, Sextans, 
et al. 1999). Low-resolution spectra of individual stars and Tj rsa Minor 
along the red giant branch (RGB) have shown that both 

Galactic and M31 satellite dSphs have internal metallic- One problem in deriving an accurate SFH from stellar 
ity dispersions (e.g., Suntzeff et al. 1993; Smecker-Hane colours is the so-called age-metallicity degeneracy. Older 
& Mc William 1999; Cote, Oke, & Cohen 1999). Using stars become redder and the same is true for more metal- 
high dispersion (R=34000) spectroscopy, Shetrone, Cote, rich stars. Hence, young and metal-rich stars show colours 
& Sargent (2001) have recently presented abundances of similar to old and metal-poor ones. We have empirically 

solved the age-metallicity degeneracy for the three dSphs 

Send offprint requests to: C. Ikuta cited above to derive their star formation histories (SFHs) 

* JSPS fellow to a level of accuracy which is impossible in other galaxies 
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in the Local Group. To reveal SFHs in detail would help 
to understand how large galaxies and the environments 
affect the evolution of dwarf galaxies. 

In addition to a study of the evolution of dwarf galaxies 
themselves, one can also learn about the possible merging 
history of the Milky Way by comparing the abundances 
of the satellite dwarf galaxies with those of Galactic field 
halo stars. In a hierarchical galaxy formation scenario, 
larger galaxies grew at the expense of smaller gaseous 
fragments (e.g., White & Rees 1978; Kauffmann, White & 
Guiderdoni 1993; Blumenthal et al. 1994; Cole et al. 1994). 
From the observational perspective, Searle & Zinn (1978) 
proposed a model in which the Galactic halo was formed 
via infall and destruction of proto-Galactic fragments. 
Increasing empirical and analytical evidence suggests that 
the Galactic halo was, at least in part, assembled from 
chemically-distinct, low-mass fragments (e.g., Searle & 
Zinn 1978; Yanny et al. 2000). Such evidence includes the 
recent discovery of the tidally disturbed dSph Sagittarius 
which is falling onto the Galaxy (Ibata, Gilmore & Irwin 
1994) as well as the numerous reports of kinematical sub- 
structure among halo field stars (e.g., Chen 1998). The sce- 
nario can be tested by comparing abundance patterns of 
field stars in the Galactic halo and those of satellite dwarf 
galaxies. If today's dwarf galaxies in the Local Group were 
the counterparts of proto-Galactic gaseous fragments, the 
abundance patterns of dwarf galaxies should be similar 
to those of field stars in the Galactic halo. The observa- 
tions by Shetrone et al. (2001), however, have shown that 
the patterns of abundance ratio in the dSphs are differ- 
ent from those of the Galactic halo stars. The dSphs have 
-0.4 < [Mg/Fe] < 0.4, while field stars in the Galactic 
halo have 0.2 < [Mg/Fe] < 0.6 . The authors concluded 
that the Galactic halo was unlikely to have formed via the 
accretion of objects similar to the dSphs Draco, Sextans, 
and Ursa Minor. 

In this paper, we discuss the SFHs and the CEHs of the 
three dSphs Draco, Sextans, and Ursa Minor. Using the 
abundance patterns and the CMDs, we show that the low 
(1—5% of the solar neighbourhood) star formation rate 
and relatively long duration (3.9 — 6.5 Gyr) of initial star 
formation explain both the abundance patterns and mor- 
phologies of CMDs such as well-populated red horizontal 
branch (RHB) and the tight red giant branch (RGB). 

The paper is organised as follows. In section 2, chemical 
evolution model is confronted with the observed stellar 
abundances to study the CEHs and SFHs. In section 3, 
by calculating CMDs by using our CMD simulation code 
(Ikuta 2001), we will show that the SFHs derived in the 
previous section succeed in reproducing observed CMDs. 
Discussions and conclusions are presented in sections 4 
and 5, respectively. 



2. Observed Abundances and Chemical Evolution 




I ■ Draco 
" aUrsa Minor 
* Sextans 



[Fe/H] 

Fig. 1. Theoretical abundance patterns together with the 
observations (Shetrone et al. 2001) for stars in the dSphs 
Draco, Sextans, and Ursa Minor. The meaning of the 
marks is written on the panel. The same IMF (the Salpeter 
IMF: x = 1.35) with upper (m u = 60M©) and lower (mi = 
0.1Af Q ) mass limits are adopted in all the model, while dif- 
ferent values of esFGyr -1 ) are assumed; model A (csf = 
2 10~ 2 Gyr" 1 : dotted line); modelB (e SF = 1 10~ 2 Gyr" 1 : 
dashed line); model C (esf = 5 10~ 3 Gyr -1 : solid line); and 
model D (esF = 2 10~ 3 Gyr -1 : dotted-dashed line). Colour 
version of this figure appears in the electric journal. 

2.1. Chemical evolution model 

Generally stellar birth rate is separated into two indepen- 
dent functions. The birth rate of stars with mass between 
to and m + dm is described as C(t)(f>(m)d'm, where C(t) 
and 0(m) are the SFR and the initial mass function (IMF), 
respectively (Tinsley 1980). The IMF is assumed to be 
time invariant with a power-law spectrum. Normalising 
to unity, we have 

where the lower and the upper mass limits are assumed 
to be 777; = 0.1M Q and 777„ = 60M Q , respectively. The 
Salpeter IMF has x = 1.35 in this definition. Similarly to 
the solar neighbourhood disc model of chemical evolution 
(Arimoto, Yoshii, & Takahara 1992), the Salpeter IMF is 
assumed in this section. Models with different IMFs will 
be discussed in section 4. The gas mass M g (t) changes 
through star formation: 

**®—o®+m <2) 

where E(t) is the gas ejection rate from dying stars, as- 
suming a closed-box model. In the case of the Galactic 
disc, an infall model is generally assumed to solve the 
so-called G-dwarf problem (e.g., Tinsley 1980; Arimoto 
et al. 1992). However, there is no observational evidence 
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Table 1. Model Parameters 
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for the G-dwarf problem in dwarf galaxies, thus we do 
not consider the infall of gas here. Gas outflow during 
the star formation is not considered, either. A recent nu- 
merical experiment showed that the fraction of mass loss 
from dwarf galaxies is small (at most ~ 7 percent of the 
galaxy mass; MacLow & Ferrara 1999). Metal-enhanced 
wind is often introduced to explain observed properties 
of dwarf galaxies (e.g., Vader 1986). MacLow & Ferrara 
(1999) showed that indeed most of the metals are able to 
leave dwarf galaxies during star burst, although the mass 
loss rate is small. Virtually no metals are retained in dwarf 
galaxies if high energy is supplied by supernova explosions 
(see table 3 in their literature for details) . Complete metal 
loss means no chemical enrichment in dwarf galaxies. This 
is inconsistent with observations which show clearly that 
dSphs are chemically polluted. Their results suggest that 
only star formation of a low rate allows a dwarf galaxy 
to retain a large fraction of heavy elements produced in 
massive stars. As will be shown later, very low SFRs are 
indeed required to explain the observed abundances and 
CMDs. Thus, a simple model is adopted here and a galaxy 
is assumed to be a closed system. 

The SFR is assumed to be proportional to the gas frac- 
tion f s (t): 

C(t) = e SF M G ■ f g (t) = e SF • M g (t) (3) 
with 

f g (t) = M g (t)/M G , (4) 

where egp (Gyr) and Mq are the timescale of star forma- 
tion and the initial total mass of a galaxy, respectively. 

The evolution of the abundance of the z-th element 
Zi{t) is given as 

= - Zl C{t) + E Zt (t) + Z t (0)M G , (5) 

where Ez, (t) is the total ejection rate of processed and 
unprocessed i-th element, Zi(0) is the initial abundance 
of i-th element in a proto-cloud. Hereafter, Z,(0) = is 
always assumed, i.e., a proto-cloud was metal- free. We cal- 
culate E(t) and Ez t (t) from the following integrals: 

r-m u 

E{t) = (1 - w m )C(t - t m )<f>{m)dm, (6) 

J m t 

EzM 

= I [(l-w m )Zi(t - t m )+pi m ]C(t - t m )(j>(m)dm, (7) 

J m t 



where t m is the lifetime of a star with mass m and the 
lower limit m t is the stellar mass with lifetime t m = t. 
Nucleosynthesis data w m and p im are the remnant mass 
fraction and the total mass fraction of processed and un- 
processed i-th. element. It is still debated whether stellar 
metallicity affects both w m and p imi but no firm consen- 
sus has been reached yet. Thus we neglect these possible 
effects: instead we adopt the values calculated for the so- 
lar abundance. We assume that the lower mass limit TOjji 
for a progenitor of Type II supernova (SNell) is 1OM 
(Tsujimoto et al. 1995). Nucleosynthesis of SNell and of 
SNela are explicitly included, since the former profoundly 
produce the a elements (e.g., O, Mg, Si, Ca), while the 
latter the iron-peak elements (e.g., Fe, Ni, Co). The data 
of nucleosynthesis are taken from Thielemann, Nomoto, 
& Hashimoto (1996) for SNell, Nomoto, Thielemann, & 
Yokoi (1996) for SNela, and Renzini & Voli (1981) for low- 
mass (< 8M©) stars. Progenitors of SNela are generally 
considered to have longer lifetimes than those of SNe II, 
since a majority of them are likely to be Chandrasekhar 
mass white dwarfs (e.g., Nomoto, Iwamoto, & Kishimoto 
1997). Observations for disc (Edvardsson ct al. 1993) and 
halo (e.g., McWilliams et al. 1995; Gratton & Sneden 
1987) stars in the Milky Way show that there is a sud- 
den decline of [a/Fe] in the [a/Fe]-[Fe/H] diagram. The 
decline is generally interpreted as an onset of SN la explo- 
sions. Studying observed [O/Fe] break at [Fe/H] ~ —1 for 
disc and halo stars in the Milky Way, Yoshii, Tsujimoto, 
& Nomoto (1996) suggested that the effective lifetime of 
SNela is around 1.5 Gyr. A new SNIa model (Hachisu, 
Kato, & Nomoto 1999) suggested that occurrence of SNela 
depends on the metallicities of progenitors as well as the 
lifetimes. Based on the new SNIa model, Kobayashi et al. 
(1998) predicted that SNIa events occur only for the pro- 
genitors with [Fe/H] > —1.1. On the observational side, 
a element to iron abundance ratios of the three dSphs, 
Draco, Sextans, and Ursa Minor (Shctrone et al. 2001) 
have indicated a significant contribution from SNela at a 
metallicity lower than [Fe/H]= —1.0. The observed trends 
of [a/Fe] are difficult to explain by the new SNIa model. 
Thus, the metallicity effect on SNela is not considered 
here and the lifetime of the progenitors is assumed to be 
1.5 Gyr (Yoshii et al. 1996). 

2.2. Comparison with the observation 

Fig. 1 reproduces the observed stellar abundances of 
the dSphs Draco, Sextans, and Ursa Minor in the 
[Mg/Fe] — [Fe/H] diagram (Shetrone et al. 2001). Despite 
the narrow widths of the RGB (see figure 1 of Shetrone 
et al. 2001), the iron abundances exhibit a large disper- 
sion (— 3 < [Fe/H] < —1.4), and the abundance ratios of 
the a-elements (e.g., Mg) to iron are near or below solar. 
The abundance ratio [Mg/Fe] starts to decrease at [Fe/H] 
~ —2. On the other hand, [Mg/Fe] in the solar neigh- 
bourhood decreases with [Fe/H] at around [Fe/H]> — 1, 
which is generally interpreted as the onset of SNela ex- 
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(1) Assumption of SFH (e SF , AT SF , Age, IMF) 






(2) Calculation of CEH 
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Fig. 2. Flow-chart of a CMD simulation code. 



plosions (e.g., Matteucci & Greggio 1986). Supposing that 
the decline of [Mg/Fc] in the dSphs is also caused by 
SNcIa, the star formation should have continued longer 
than the lifetime (typically 1 — 2Gyr) of progenitors of 
SNela and the SFRs should be much lower than that of 
the solar neighbourhood. A time-delay model (Matteucci 
& Brocato 1990) of the SNela precisely predicted that the 
low star formation rate results in lower [Mg/Fe] ratios rel- 
ative to stars in the Galaxy. 



In Fig. 1, theoretical abundance patterns of our models 
A— D are superposed. The models assume different £sf of 
2 • 1CT 2 , 1 • 1CT 2 , 2 • 1CT 3 , 1 • 1(T 3 (Gyr' 1 ) as indicated in 
Fig. 1. The model parameters are summarised in table 1. 
A typical value of egp for the solar neighbourhood disc is 
e S F = 0.2 Gyr" 1 (Arimoto et al. 1992). 

The models B— D shown in Fig. 1 give a good fit to the 
observed abundance patterns. Since the abundance ratios 
resulting from the model A are too high to reproduce the 
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Fig. 3. CMD of the CO region of the Draco dSph observed 
by Piatek et al. (2001). The points marked with open cir- 
cles correspond to some confirmed RR Lyrae stars (Baade 
& Swope 1961). 

observed ones, the model A is rejected. The model D with 
csf = 1 • 10~ 3 Gyr -1 is not appropriate, either. Since 
chemical evolution is so slow, it takes around 15 Gyr to 
reach [Fc/H]= —1.4, which is the highest abundance ob- 
served in these galaxies. If this model is adopted, the star 
formation must continue till now. Clearly this is inconsis- 
tent with the fact that no star formation occurs currently 
in these galaxies. Thus, only models B and C remain as 
being acceptable. We note that the observational data for 
the dSph Sextans have lower signal-to-noise ratios than 
those for the dSphs Draco and Ursa Minor (see table 2 
of Shetrone et al. 2001) and the abundances of the dSph 
Sextans are less reliable than those of the other two dSphs. 
Thus, it is not necessary to take the discrepancy between 
the model predictions and the abundances of the dSph 
Sextans too seriously. 

In short, the abundance patterns of the dSphs Draco, 
Sextans, and Ursa Minor all suggest that the chemical en- 
richment occurred with SFRs which are much lower than 
normal spiral galaxies (only 1 — 5%). 

3. CMDs 

The models B and C fit well to the observed trends in the 
[Mg/Fe] — [Fe/H] diagram. Next, we discuss if the colour- 
magnitude diagrams (CMDs) predicted by these models 
can reproduce the observed ones. The CMDs provide us 
with the most detailed information to follow SFHs and 
CEHs back to the oldest stars. The SFH derived from the 
abundance patterns in the previous section must give a 
consistent CMD with the observations. 



3.1. Simulation code of a CMD 

Simulating numerical CMDs has become a standard tech- 
nique to study the SFHs of nearby galaxies through the 
observed CMDs. Galaxies are composed of complex stellar 
populations, and the observed CMDs are affected by pho- 
tometric errors and detection incompleteness. A Monte 
Carlo simulation allows us to take into account both si- 
multaneously. A composite stellar population is randomly 
generated according to predictions of stellar evolution the- 
ory with the IMF and SFR being assumed, and then un- 
certainties of the data such as the increasing scatters and 
rising incompleteness at fainter magnitudes are taken into 
account. 

Because of these advantages, this approach is becom- 
ing more and more frequently used to study the stellar 
populations of nearby galaxies. The approach was first 
applied by Ferraro et al. (1989) and was fully described 
by Tosi (1991) and Greggio et al. (1993). A more quanti- 
tative method, the R-method, to compare simulated and 
observed CMDs was presented by Bertelli et al. (1992) 
and was described by Vallcnari et al. (1996) in more 
detail. The statistical comparison was further developed 
by Tolstoy & Saha (1996) and Tolstoy (1996) by adopt- 
ing Bayes' theorem. Now several groups have constructed 
CMD simulators and have investigated SFHs of galaxies 
in the Local Group (e.g., Aparicio et al. 1996; Dolphin 
1997; Hernandez, Valls-Gabaud, & Gilmore 1999; Gallart 
et al. 1999). 

We should point out, however, that there were two 
problems in the previous studies. First, all the approaches 
so far developed use the so-called optimising method of 
stellar population synthesis. The best mixture of stellar 
population is searched iteratively to reproduce the CMDs. 
Problems with this approach are the unproven uniqueness 
of the solution (e.g., Greggio et al. 1999) and a lack of 
evolutionary information. 

Second, stellar metallicity was assumed a priori in 
all the simulation codes presented previously. Some in- 
troduced metallicity variation in time (e.g., Gallart et al. 
1999), but the variation assumed was independent of the 
SFH. Since the metallicity affects evolutionary tracks and 
atmospheres of stars, the colours and luminosities of stars 
are changed by metallicity. Thus, it is crucial to use a 
simulator of CMDs which fully takes into account a CEH 
to interpret the CMDs properly and to derive the SFHs 
accurately. The simplification and/or neglect of chemi- 
cal evolution lead to a serious problem called the age- 
metallicity degeneracy. Since stellar colours become bluer 
when stars are younger and/or poorer in metallicity, there 
are at least two interpretations for a given position of a 
star in the CMD; young and metal- rich or old and metal- 
poor. Therefore, an assumption of the metallicity (or age) 
of a certain stellar population may result in wrong esti- 
mation of age or metallicity, due to the age-metallicity 
degeneracy. 

Aiming to solve the degeneracy, we have built a numer- 
ical simulation code of CMD morphology (Ikuta 2001). In 
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Fig. 4. Simulated CMDs of models B (left) and C (right). In modelB, e SF = 1 1(T 2 (Gyr" 1 and AT = 3.9 Gyr are 
assumed, while egp = 5 1CP 3 (Gyr -1 and AT = 6.5 Gyr are adopted in model C. In the both models, we assumed the 
same upper and lower mass limits (m u = 60M Q and mi — O.1M , respectively), the Salpeter IMF, and the galactic 
age of 12 Gyr. Observational errors and data incompleteness are taken from Piatek et al. (2001). 



the code, we adopt an evolutionary method of stellar pop- 
ulation synthesis (e.g., Arimoto & Yoshii 1986). A galaxy 
is assumed to have been a proto-galactic gas cloud at the 
beginning. Stars formed and newly processed elements 
were ejected from stars at the end of their lives either 
via stellar wind or via supernova explosions (SNela and 
SNell). The gas was chemically enriched and the next gen- 
eration of stars changed their evolution and photometric 
properties due to a metallicity increase. The code is partic- 
ularly designed to disentangle stellar age and metallicity 
both of which heavily affect the morphology of CMDs. 
To take into account effects of chemical evolution, a fine 
interpolation in stellar ages and metallicities is adopted 
in our code. The later stages of stellar evolution, in par- 
ticular, horizontal branch (HB), asymptotic giant branch 
(AGB), and post-AGB stars are treated in a sophisticated 
way, since these bright stars are crucial for this purpose. In 
this section, the procedure of our Monte Carlo simulation 
of CMDs and the ingredients are detailed. 

3.1.1. Outline of the Simulation Procedure 

Fig. 2 shows the flow-chart of the simulation. First, a SFH 
is assumed. The parameters to describe the SFH are the 
timescale of star formation e<7p (Gyr), the duration of star 
formation ATsf (Gyr), the slope of the IMF (x), and the 
age of a galaxy. Second, chemical evolution is calculated, 
which gives the total number of stars at each metallicity at 
each time step. Third, a HR diagram is calculated by us- 
ing stellar evolutionary tracks, where a Monte Carlo sim- 
ulation is performed to assign stellar masses. Fourth, the 



effective temperature T c g and luminosity L of a star are 
transformed to colour and magnitude by using a library of 
stellar model atmospheres. Fifth, a Monte Carlo simula- 
tion is performed again to take into account observational 
errors and detection completeness, and thus the CMD is 
obtained. A number of stars in the simulated CMD is ad- 
justed to an observed CMD. Finally, the simulated CMD 
is confronted with the observed one by visual inspection. 
If they are inconsistent, the parameters of the SFH are 
changed, and the simulation is re-started. If the simulated 
CMD fits the observation, the SFH and CEH are derived. 
Ingredients of the simulation, such as stellar evolutionary 
tracks, model atmospheres, and observational conditions, 
are described in the subsequent section. 

3.1.2. Stellar Tracks and Model Stellar Atmosphere 

Stellar evolutionary tracks give luminosities, effective tem- 
peratures, and surface gravities of stars with given mass 
and chemical composition as a function of age. The tracks 
depend on the basic parameters such as the initial mass 
and chemical compositions. The properties of stellar popu- 
lations in galaxies also depend on other properties not ex- 
plicitly included in most of current stellar evolution mod- 
els, i.e., the stellar rotation and close binary companions. 
The numerical calculations are also sensitive to the treat- 
ment of the convection, such as a mixing length parameter 
I /H p and convective overshooting. 

The Padova stellar evolutionary tracks are adopted for 
stars from the main sequence (MS) to the early asymptotic 
giant branch (EAGB) stars. The Padova tracks cover wide 
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Fig. 5. Simulated CMDs of models B (panel (a)) and C (panel (b)). Adopted errors and data completeness are the same 
as those in Fig. 4, while the number of stars in each panel is increased for presentation purposes. Colours correspond 
to the following metallicity ranges; blue: [M/H]< -2.0; green: -2.0 < [M/H] < -1.8; orange: = 1.8 < [M/H] < -1.6; 
red:-1.6 < [M/H]. 



ranges in age (0 to 16 Gyr) and metallicity (Z = 0.0001 to 
Z = 0.05), and are one of the most complete sets of stel- 
lar evolutionary models currently available. The Padova 
tracks were calculated with revised radiative opacities 
(Iglesias, Rogers, & Wilson 1992) and with l/H p — 1.63. 
The adopted tracks were calculated for the following set of 
chemical compositions: (Y — 0.230, Z = 0.0001; Girardi 
et al. 1996), (Y = 0.230, Z = 0.0004; Fagotto et al. 1994a), 
(Y = 0.240, Z = 0.004; Fagotto et al. 1994b), (Y = 
0.250, Z = 0.008; Fagotto et al. 1994b), (Y = 0.280, Z = 
0.02; Bressan et al. 1993), and (Y = 0.352, Z = 0.05; 
Fagotto 1994a). 

To obtain the isochrone for a given metallicity, the 
original isochrones are linearly interpolated in metal- 
licity. Since the metallicity changes stellar evolution 
non-linearly, any extrapolation in metallicity should be 
avoided. Thus, stars of metallicity lower (Z < 0.0001) and 
higher (Z > 0.05) than the Padova tracks do not appear 
in our simulated CMDs and are not considered here. This 
is justified, because no significant populations of such ex- 
treme metallicities are known to exist in dwarf galaxies. In 
addition to these, close binary systems are not included, 
since the input tracks are only for isolated single stars. 
Although our code does not explicitly include binary stars, 
it effectively takes into account apparent binary stars as 
stellar blends. Gallart et al. (1999) reported a study that 
is more detailed on this issue. 

For the RGB evolution, the mass loss law of Reimers 
(1977) is adopted; 



M = 



-4- 10~ 13 r? 



L/L & 



(MQ-yr- 1 ), 



or equivalently, 

log(-M) = - 4.67 + log?; + 1.51og(L/L ) 

-log(M/M Q ) - 21ogT cff , (9) 

where L, g, R, T c ff , and M are stellar luminosity, surface 
gravity, radius, effective temperature, and initial mass, re- 
spectively. The mass loss efficiency parameter 77 is assumed 
to be 1/3 (Renzini & Voli 1981). 

The horizontal branch (HB) stars are distributed on 
the HR diagram according to a modified Gaussian mass 
distribution equation (Lee et al. 1990); 



$(M) cx [M - (M HB - AM)] (M RGB - M) 



x exp 



(M HB - Mf 



(10) 



(8) 



where a is a mass dispersion factor, and Mhb = (-MrGB — 
AM) is the mean mass of HB stars. A value of u = 
O.O6M has been chosen to represent the HB distribu- 
tions similar to those of Galactic globular clusters (Rood 
1973; Lee et al. 1990). This is equivalent to considering the 
statistical dispersions of mass loss rate along the RGB. 

To convert the theoretical temperature-luminosity 
data to the observable colour-magnitude plane, a stellar 
spectral library by Lejuene, Cuisinier & Buser (1998) is 
used. This consists of Kurucz's (1992) spectra for hot- 
ter stars (O-K), Bessell et al.'s (1989; 1991) and Fruks et 
al's (1994) spectra for M giants, and Allard & Hauschudt's 
(1995) for M dwarfs. In the original model spectra, system- 
atic deviations become apparent when colour-temperature 
relations computed from the models are compared to the 
empirical ones at IZq. The library adopted here is a 



8 



C. Ikuta & N. Arimoto: Extended Star Formation in Satellite Dwarf Galaxies 



version that the authors made by correcting the origi- 
nal library. The corrections are especially important for 
M star models. The fundamental parameters (T e g, g, and 
[M/H]= log(Z/Z Q )) are wide enough to cover all spec- 
tral types and luminosity classes that appear in observed 
CMDs; 3 x 10 3 < T cff < 5 x 10 4 , -1.02 < \ogg < 5.5; 
—5.0 < [M/H] < 1.0. To reduce computational time, syn- 
thetic magnitudes and colours tabulated by Lejuene et al. 
(1998) are linearly interpolated in logZ, logT c ff , and logg. 

3.2. Comparison with the CMD of the Draco dSph 

In section 2, models B and C are shown to fit well the 
observed trends in the [Mg/Fe] — [Fe/H] diagram. Next, 
we discuss if the CMDs predicted by these models can 
reproduce the observed ones. Ideally, the SFH would be 
studied using a CMD which satisfies the following two 
conditions: (1) it covers a whole galaxy; (2) it is signifi- 
cantly deeper than the turnoff of the oldest stellar popula- 
tion. Unfortunately, such CMDs are not available for the 
galaxies discussed here. The situation, however, is better 
for the dSph Draco. Recently Piatek et al. (2001) imaged 
nine fields in and around the Draco dSph by using the 
KPNO 0.9m telescope and presented the CMDs down to 
a luminosity level ~ 2 mag fainter than the HB. Thus, we 
mainly discuss the dSph Draco in this section and briefly 
mention the dSphs Sextans and Ursa Minor. 

Fig. 3 shows the observed CMD of the CO field in the 
Draco dSph (Piatek et al. 2001), which is characterised 
by the narrow RGB and the well-populated red HB. The 
HB in figure 1 of Piatek et al. (2001) is somewhat dis- 
torted due to RR Lyrae stars which are marked with open 
circles in Fig. 3. Therefore, we do not take seriously the in- 
consistency between the observed and simulated blue HB 
morphology (V— R< 0.1). Historically, Baade & Swope 
(1961) obtained the CMD of the Draco dSph and found 
that the RGB is generally similar to those of metal-poor 
globular clusters, although it is rather wide. The authors 
also found that the populous RHB is incompatible with 
the low metallicity of a Population II system. This is the 
notorious second-parameter problem of the dSph Draco. 

Figs. 4a and 4b represent the CMDs simulated by mod- 
els B and C, respectively, where the galactic age is as- 
sumed to be 12 Gyr, which is similar ages of Galactic 
globular clusters (e.g., Carretta et al. 2000); values of other 
parameters are written in table 1. The photometric errors 
and detection completeness in the simulations are taken 
from Piatek et al. (2001). Both models can reproduce char- 
acteristics of the CMD of the dSph Draco, such as narrow 
RGB and heavily populated RHB. Despite the metallicity 
dispersion (-3 < [Fe/H] < -1.4), the theoretical CMDs 
results in narrow RGBs. This is due to the so-called age- 
metallicity degeneracy as explained below. 

Fig. 5 shows CMDs simulated by assuming models B 
and C, where different colours represent different metal- 
licity range as described in the figure. For presentation 
purpose, Figs. 5a and 5b contain larger numbers of stars 
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Fig. 6. The same as Fig. 1, but for models E-G. Different 
lines mean different IMF powers; model E (x — 1.75: 
dashed line); modclF (x = 1.95: solid line); modelG 
(x — 2.15; dotted-dashed line). The star formation time 
scale is the same as for the solar neighbourhood disc, i.e., 
£sf = 0.2 (Gyr -1 ). See electric edition of the journal for a 
colour version of this figure. 



than those of Figs. 4a and 4b, respectively. The adopted 
observational conditions are the same as those of Fig. 4. 
Clearly, there is little correlation between the colour and 
metallicity of RGB stars because of the age-metallicity 
degeneracy. Since stars of higher metallicity are younger, 
the metallicity and age effects on the RGB colours are 
cancelled out. This keeps the RGB narrow and tight. The 
age effect appears in the HB morphology too. Since only 
metal-poor and old (> 10 Gyr) stars can evolve to BHB, 
younger core-He burning stars lie in the red part of the 
HB. The populous red HBs in Fig. 5 clearly show this age 
effect on the HB morphology. A weak correlation between 
the stellar colours and metallicities is found in the bright 
(m v < 18 or M v < -1) RGB. 

In model C, the RHB is more populous and brighter 
than in model B. This is because the longer duration of 
the star formation period increases the number of younger 
and more massive core-He burning stars. However, the 
difference between models B and C is so small that it 
is difficult to determine the final epoch of star formation 
from the CMD that does not reach down to the turnoff 
level. Nevertheless, it can be said that, in the dSph Draco, 
chemical evolution was very slow and the duration of star 
formation was at least longer than 3.9—6.5 Gyr. 

This evolutionary picture can explain both the narrow 
RGB and the populous red HB of the Draco dSph. Thus, 
the second parameter problem in the dSph Draco is solved 
for the first time by introducing relatively long duration of 
star formation (> 3.9 — 6.5 Gyr) which is fully consistent 
with the observed abundance patterns. 
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Fig. 7. The same as Fig. 4, but for models F and G. 



4. Discussion 

The abundance patterns and the CMDs of the dSphs 
Draco, Sextans, and Ursa Minor share the common fea- 
tures. It can be said that they have similar SFHs and 
CEHs. The CMD of the Sextans dSph (see figure 2 in 
Suntzeff et al. 1993; figure 1 in Shetrone et al. 2001) is 
characterised by a narrow RGB and a red HB similar to 
the dSph Draco. This suggests that the SFH of the dSph 
Sextans is similar to that of the dSph Draco. In the CMD 
of the Ursa Minor dSph (see figure 1 in Shetorne et al. 
2001), the blue HB is more populous and the RGB is nar- 
rower. This might suggest that the time scale of star for- 
mation was shorter in the dSph Ursa Minor. Nevertheless, 
it can be said that, in the Draco dSph, the SFR was very 
low and the duration of star formation was at least longer 
than 3.9—6.5 Gyr. For the other dSphs, it is safe to say 
that the characteristics of the SFHs are the low SFRs and 
the long durations (> several Gyrs) of the star formation 
period. 

For the dSph Draco, model C may conflict with the 
CMD obtained by the Hubble Space Telescope (Grillmair 
et al. 1998). Through a comparison with the fiducial lines 
of metal-poor globular clusters M68 and M92, Grillmair 
et al. (1998) suggested that the dSph Draco is older than 
M68 and M92 by 1.6 ± 2.5 Gyr. Their result, however, 
seems inconclusive. Firstly, the sparseness of the photo- 
metric sample makes it difficult to measure the contribu- 
tion of any intermediate-age population. Although a pre- 
dominant RHB is the characteristic of the dSph Draco, the 
few RHB stars appeared in their CMDs. Secondly, they 
used different techniques for photometry and presented 
the combined results. For bright stars aperture photom- 
etry was adopted, while the point spread function fitting 
method was used for faint stars. The border is unclear. 



This causes systematic discrepancies between the photom- 
etry of faint and bright stars. Third, a comparison with 
M68 and M92 is questionable. They used the F606W filter 
of the WFPC2 system. However, the WFPC2 Instrument 
Handbook recommends that the F555W filter is a better 
approximation to the Johnson V-band than the F606W. 
Their filter selection could have reduced the accuracy of 
the photometric calibration. Because of the combination 
of these uncertainties, it is dangerous to place much faith 
on their age estimation. More data are required to clarify 
these matters and images of a wide special coverage are 
particularly important. 

In the above discussion, the decline in [Mg/Fe] ob- 
served in the dSphs is interpreted as the effect of the 
chemical enrichment by SNela. Since the signal-to-noise 
ratios of the data (Shetrone et al. 2001; S/N~ 13 - 36) 
for the dSphs are lower than other observations for metal- 
poor stars (e.g., McWilliams et al. 1995;S/N~ 30 - 40), 
one might argue that the quality of the data is not high 
enough to study the CEHs. However, the goal of Shetrone 
et al. (2001) was the measurement of overall abundance 
differences, not absolute values as presented in their study, 
and it can be said that [a/Fe] declines at [Fe/H]~ —2 on 
average (see figures 4 and 5 in Shetrone et al. 2001). Thus, 
it is a pertinent interpretation that enrichment by SNela 
causes the decline of [a/Fe]. Spectroscopic observations of 
high signal-to-noise ratios arc required to confirm this. 

We have shown that low SFRs can explain the observed 
trends in the [Mg/Fe] — [Fe/H] diagram based on models 
assuming the Salpeter IMF. However, models of steeper 
IMFs and higher esf can also reproduce the trends. The 
steeper IMFs lead to slow chemical evolution, so that 
SNela start to explode at low metallicity. This results in 
the decline of [Mg/Fe]. 
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Fig. 8. Age-metallicity relations of models A— D (panel (a)) and E— G (panel (b)). The adopted values of csf (Gyr" 1 ) 
or IMF power are shown on the panel. Colour version of this figure appears in the electric edition of the journal. 



Fig. 6 shows chemical evolution models E— G assuming 
different IMFs, where the IMF slopes are x = 1.75, 1.95, 
and 2.15, respectively, and £sf = 0.2 (a standard SFR 
for the solar neighbourhood; Arimoto et al. 1992). These 
models, except for model E, are also consistent with the 
observations. As with model A, the abundance ratios pre- 
dicted by model E are too high and thus inconsistent 
with the observations. For the observed metallicity range 
([Fe/H]< -1.4), the trends of [Mg/Fe] of models E-G 
are very similar to those of models A— D, although they 
become different at higher metallicity (see Fig. 1). 

Figs. 7a and 7b represent the CMDs simulated by 
models F and G, respectively. The photometric errors and 
detection completeness are taken from Piatek et al. (2001), 
i.e., the same as the simulations of models B and C shown 
in Fig. 4. The CMDs in Fig. 7 resemble those in figure 4 
in terms of the morphologies of RGB and HB. The slow 
chemical evolution and the long (> 1.6 Gyr) durations 
of star formation of models F and G result in a narrow 
RGB and the populous RHB by the same reasons as those 
of models B and C. This indicates that both abundance 
patterns and morphologies of the RGB and HB become 
similar when models have similar age-metallicity relations 
(AMRs). 

Figs. 8a and 8b show the AMRs of models A-D 
and E— G, respectively. The figures demonstrate that the 
AMRs of models E—G with the steeper IMFs result in very 
similar to those of models A— D with the lower SFRs at 
low- metallicity ([Fe/H]< —1.4). Since both the abundance 
patterns and the morphologies of the RGB and HB are af- 
fected by CEHs, an additional information is necessary to 
solve the degeneracy of IMF and SFR. In the models E— G 
the initial star formation lasts only 1.6 — 2.2Gyrs, while 
in the models B— C it continues as long as 3.9 — 6.5 Gyrs. 
Therefore, the models B— C should give systematically 
brighter turnoff magnitudes than the models E—G, and 
magnitude difference between the main sequence turnoff 
and the HB F(TO - HB) in the models E-G should be 



smaller than in the models B—C. The models E—G give 
^(TO - HB) is - 3.1 mag, while the models B and C 
does V(TO — HB) ~ 3.1 mag and ~ 2.8 mag, respectively. 
Deeper images with wider field of views will allow us to 
determine the final epochs of the initial star formation. 

Figs. 1 and 6 show that [Mg/Fe] converges at higher 
metallicity in the models with the different SFRs and 
that it does not in those with the different IMFs. This 
difference allows us to derive the IMF and SFR indepen- 
dently. Abundance patterns in more massive dwarf galax- 
ies and/or dwarf irregular galaxies will test whether the 
IMF or SFR produces in the abundance patterns and 
low metallicities in the dSphs. We note that even steeper 
(x > 2.75) IMFs are required to explain the low metallic- 
ities of dwarf irregular galaxies if the IMF is the primary 
cause of the low metallicities. The CMD of the dwarf irreg- 
ular galaxies in the Local Group reveal that they contain 
very old (> 10 Gyr) populations as well as young ones, i.e., 
the star formation and chemical enrichment have contin- 
ued for at least 10 Gyr. Since they are still metal-poor, 
extremely steep (x > 2.75) IMFs are required, which have 
never been reported. Therefore, we believe that low SFRs 
are the primary characteristic of the SFHs of the dSphs. 

So far, the gas infall and outflow during star forma- 
tion have been neglected, since no evidence of the gas 
infall and outflow has presented. Metallicity distribution 
could provide a clue to judge whether gas infall or out- 
flow should be considered. We stress that spectroscopic 
observations is indispensable to derive a metallicity dis- 
tribution. Fig. 5 demonstrates little correlation between 
the colours and metallicities of the RGB stars. Because 
of the age-metallicity degeneracy and contamination by 
AGB stars, a stellar colour is not always a good indica- 
tor of the metallicity. Thus, a careful procedure should be 
adopted to convert stellar colours into metallicities (Ikuta 
2001), although a relatively simple technique is often used 
(e.g., Harris & Harris 1999). 
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Long lasting star formation at a very low rate explains 
the observed trend of [Mg/Fe] and the CMDs. Because of 
the low SFRs, plenty of gas (~ 97 percent of the galaxy 
mass) still remain even at the final epoch of star formation 
(i.e., t — ATsf). The remaining gas has to be removed 
to complete star formation and to evolve to a gas-poor 
system. Energy injection from supernovae could not be a 
sufficient mechanism to expel the gas from dwarf galaxies 
(MacLow & Ferrara 1999; Ikuta 2001). If this is the case, 
gas removal from dwarf galaxies should result from ex- 
ternal mechanisms such as ram pressure stripping and/or 
tidal shocks. Studying the SFHs of 18 dwarf galaxies in 
the Local Group based on the CMDs derived by a uni- 
form method of photometry, Ikuta (2001) found that the 
durations of star formation period correlate with the dis- 
tances from the Milky Way or M31. The correlation may 
suggest that the environmental effects play a key role in 
the evolution of the dwarf galaxies. Numerical simulations 
(Moore et al. 1998; Mayer et al. 2001) clearly showed that 
late-type dwarf galaxies entering the dark matter halo of 
a massive galaxy are transformed into early-type owing to 
repeated tidal stripping and dynamical instabilities. Thus, 
we conclude that the star formation in the dSphs Draco, 
Sextans, and Ursa Minor terminated due to gas stripping 
by the Milky Way. 

The hierarchical clustering galaxy formation model 
suggested that star formation in dwarf galaxies can only 
occur before the cosmological re-ionisation epoch and 
there is no major star formation activities later because 
gas cannot cool and condense to form stars (e.g., Cen 
2001). This appears to be inconsistent with stellar popula- 
tions observed in the Local Group dwarf galaxies. CMDs 
obtained by the HST (e.g., Ikuta 2001) and our study 
presented in this paper clearly show extended and recent 
star formation activities in the Local Group dwarf galax- 
ies. The issue is still being debated. However, a more recent 
simulation (Kitayama et al. 2001) suggested that star for- 
mation can occur in small objects if they have baryonic 
mass larger than the threshold mass of 10 9 Af Q and 1O 6 M 
at redshifts of < 3 and ~ 5, respectively. The masses of 
today's dSphs in the Local Group are 10 7 — 10 8 M Q (e.g., 
Mateo 1998) and our results imply that progenitors of to- 
day's dSphs lost ~ 97 % of their mass after the long lasting 
star formation. Therefore, the masses of the progenitors 
are estimated at 10 9 — 10 10 M Q . Since these exceed the 
threshold predicted by Kitayama et al. (2001), the SFHs 
derived here do not contradict recent simulations of the 
formation of dwarf galaxies. 

5. Conclusion 

Based on comparisons between the theoretical chemical 
evolution models and the observed abundance patterns, 
we conclude that the initial star formation continued for 
a long duration (> 3.9 — 6.5 Gyr) in the dSphs Draco, 
Sextans, and Ursa Minor. Our simulation of the CMDs 
shows that the long duration of star formation can solve 
the second parameter problem of the Draco dSph. Because 



of the age-metallicity degeneracy, the RGBs are kept nar- 
row and tight despite their large metallicity dispersions. 

We have discussed the two cases which are consistent 
with the observed abundance patterns and the CMDs. The 
first case is a combination of low SFRs (egp = 1 • 10~ 2 — 5 • 
10" 3 Gyr" 1 ) and the Salpeter IMF (x = 1.35), while the 
second is a combination of the solar neighbourhood SFR 
(e SF = 0.2 Gyr" 1 ) and the steeper IMFs (a; = 1.75-2.15). 
The two cases are discriminated neither by the abundance 
patterns nor by the CMDs of bright stars. However, the 
initial star formation period is as long as 3.9 — 6.5 Gyrs if 
the SFR is low, while it is much shorter (1.6 — 2.2 Gyrs) 
if the SFR is high. Thus, turnoff magnitude should be 
different between the two cases. Deeper images with wider 
field of views allow us to determine the final epochs of the 
initial star formation. 
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